Indication for Primordial Anisotropics in the Neutrino Background from WMAP 

and SDSS 
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We demonstrate that combining Cosmic Microwave Background anisotropy measurements from 
the 1st year WMAP observations with clustering data from the SLOAN galaxy redshift survey yields 
an indication for primordial anisotropies in the cosmological Neutrino Background. 



Introduction Recent cosmological data coming 
from measurements of the Cosmic Microwave Back- 
ground (CMB) anisotropies (see e.g. 0), on galaxy 
clustering (see e.g. "5]) and, more recently, on Lyman- 
alpha Forest clouds (see e.g. 0) are in spectacular 
agreement with the expectations of the so-called stan- 
dard model of structure formation, based on primor- 
dial, purely adiabatic inflationarjM^erturbations and 
a cosmological constant (see e.g. [3,0). 

Cosmology is therefore becoming more and more a 
powerful laboratory where physics not easily accessi- 
ble on Earth can be tested and verified. An excel- 
lent example of this comes from the new cosmologi- 
cal constraints on neutrino physics. According to the 
standard model, a cosmic Neutrino Background (NB), 
similar to the CMB, should exist in our universe with 
a density of approximately n^, 57/cm'^ per flavor. 
Neutrinos of mass ^ 10~^eV would still be relativis- 
tic today with a Fermi-Dirac spectrum at a tempera- 
ture T ~ 2K, while more massive neutrinos would be 
non relativistic and clustered around galaxies with a 
typical velocity of ~ 200km/s. Therefore, despite 
the high density of cosmological neutrinos a direct de- 
tection is virtually impossible due to their extremely 
low energy and cross sections (see e.g. 0). 

However, cosmological neutrinos have a profound 
impact on cosmology since they change the expansion 
history of the universe and affect the growth of per- 
turbations (see jg for a detailed account). As a con- 
sequence, recent cosmological data do provide strong 
- albeit indirect - evidence for the presence of a NB 
(see e.g. and have been used to put upper limits 
on absolute neutrino masses competitive with those 
from laboratory experiments (see e.g. 0). 

In this letter we show that current cosmological data 
now provide for the first time an interesting indication 
for primordial anisotropies in the NB. Although infla- 
tionary anisotropies in the NB at the level of ~ 10~^ 
are expected in the standard scenario, a direct detec- 
tion is clearly impossible. However, anisotropies in 
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FIG. 1: The effect of NB anisotropies on the CMB temper- 
ature angular power spectrum. The standard model with 
parameters which provide the maximum likelihood best fit 
to WMAP is plotted (see jj) against the same model but 
with no NB anisotropies. The 1st year WMAP data is also 
plotted for comparison. 



the NB background affect the CMB anisotropy an- 
gular power spectrum at level of ~ 20% through the 
gravitational feedback of their free streaming damping 
and anisotropic stress contributions 9] and an indirect 
detection is indeed possible. 

Data analysis A way to parameterize the 
anisotropies in the NB has been introduced in 
with the "viscosity parameter" c^j^, which controls 
the relationship between velocity/metric shear and 
anisotropic stress in the NB. A value of c^j^ = 1/3 
is what one expects in the standard scenario, where 
anisotropies are present in the NB and approximate 
the radiative viscosity of real neutrinos (see Fig. 1). 
The case cjjg = 0, on the contrary, cuts the Boltzmann 
hierarchy of NB perturbations at the quadrupolc, forc- 
ing a perfect fluid solution with no NB anisotropies 
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but only density and velocity (pressure) perturba- 
tions. Furthermore, the generation of polarization is 
also modified, due to changes induced in the polar- 
ization source term. Cosmological data which would 
discriminate c^j^ = in favor of c^-^ — 1/3 would 
therefore provide a strong indication for the existence 
of NB anisotropics, as argued in 11]. 

We show that this is indeed the case by computing 
the constraints on c^jg obtained by combining several 
cosmological data sets. The method we adopt is based 
on the publicly available Markov Chain Monte Carlo 
package cosmomc^ |12|. which has been modified to al- 
low for values of c^jg / 1/3. We sample the following 
7 dimensional set of cosmological parameters, adopt- 
ing flat priors on them: c^j^, the physical baryon and 
CDM densities, ojb = ^bh^ and ojc = ^ch^, the ratio 
of the sound horizon to the angular diameter distance 
at decoupling, 0^, the scalar spectral index and the 
overall normalization of the spectrum, Us and Ag, and, 
finally, the optical depth to reionization, r^. Further- 
more, we consider purely adiabatic initial conditions, 
we impose flatness and we do not include gravitational 
waves. We first restrict our analysis to the case of 3 
massless neutrino families, thereby fixing the back- 
ground neutrino density, but we proceed to relax this 
assumption at the end. Introducing a neutrino mass 
in agreement with current neutrino oscillation data 
does not change our results in a significant way. 

We include the first-year data T^l (temperature and 
polarization) with the routine for computing the like- 
lihood supplied by the WMAP tear n Il4| . as well as 
the CBI da], VSA [3 and ACBAR [13^3 measure- 
ments of the CMB. The MC conver genc e diagnostics 
is done along the lines described in In addition 
to the CMB data, we also consider the constraints on 
the real-space power spectrum of galaxies from the 
SLOAN galaxy redshift survey (SDSS) '20I. We re- 
strict the analysis to a range of scales over which the 
fluctuations are assumed to be in the linear regime 
{k < 0.2ft.~"'^Mpc). When combining the matter power 
spectrum with CMB data, we marginalize over a bias b 
considered as an additional nuisance parameter. Fur- 
thermore, we make use of the HST measurement of the 
Hubble parameter Hq = lOOh km s^^Mpc"^ 21] by 
multiplying the likelihood by a Gaussian likelihood 
function centered around h — 0.72 and with a stan- 
dard deviation cr = 0.08. 

In Fig. 2 we plot the constraints obtained from our 
analysis in the c^-^^ — Us and c^^^ — Ag planes. Including 
CMB and matter power spectrum data, the best fit for 
the standard case with c^jg = 1/3 has a = 1482.9, 
while for the cjjg = case we obtain = 1490.2, 
clearly favoring the standard case (Ax^ = 7.3). Fur- 
thermore, the best fit model with c^jg = has a rather 
tilted spectrum (n^ = 0.90), while the combined effect 
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FIG. 2: Joint 2-dimensional posterior probability con- 
tour plots in the 0^;^ — (top) and c^jg — As (bot- 
tom) planes, showing the 68% and 95% contours from 
the WMAP+other CMB data alone (cyan/light gray) 
and adding SLOAN matter power spectrum information 
(black). Degeneracies between these parameters are evi- 
dent and no lower limit on c^jj, can be placed with CMB 
data alone. 



of a red index and a slightly smaller amplitude of the 
fluctuations reduces the optical depth to — 0.07. As 
already noticed in [TJ, it is in fact possible to iden- 
tify a correlation between the amplitude As and spec- 
tral index ris of the primordial fluctuations and the 
value of c^ig. This arises because the absence of NB 
anisotropics boosts the amplitude of the CMB acous- 
tic peaks, see Fig. 1. This can be compensated, for 
example, with a lower value for the spectral index n^. 
The degeneracy is nearly exact for CMB anisotropics 
alone and we found that no strong constraint can be 
placed from the CMB data only: in that case, the best 
fit value for the standard case is — 1452.5, while 
the c^ig = case has — 1454.5. However, the latter 
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FIG. 3: Illustration of the role of galaxy clustering data 
in constraining c^jg. With CMB data only (scatter plot), 
lower values of cjjg can be compensated by a redder spec- 
tral index and a lower amplitude of the fluctuations. In- 
clusion of the SLOAN data disfavors models with low Us, 
thereby cutting away most samples with low c^jj,. The 
solid black contour is the joint 2a posterior from CMB 
and SLOAN. 



case has as best fit values h = 0.79 and = 0.89, 
which are somewhat at odds with the HST measure- 
ment and with the combined WMAP-t-SLOAN anal- 
ysis, which indicates a spectral index close to scale 
invariance. The value of c^jg has very little impact 
on the matter power spectrum, but inclusion of the 
SLOAN data does reduce the allowed range of the 
other cosmological parameters, and especially of the 
tilt, thereby allowing to place more stringent bounds 
on c^ig, as illustrated in Fig. 3. The ID marginal- 
ized posterior distribution for c^jg yields a lower limit 
c^ig > 0.12 (2cr C.I., 1 tail). The value c^j^ = is 
found to lie about 2Aa away from the standard case, 
c^ig — 1/3. Therefore we can conclude that the case 
where the NB does not have anisotropics above the 
first moment is quite clearly disfavored. 

Finally, it is interesting to relax the assumption 
on the number of massless neutrino species, N^, in 
order to test the prediction of the Standard Model 
that (c^jg, A^^) ~ (1/3,3). Departures from this point 
would constitute a strong indication for the presence 
of new physics. By treating N^, as a free parameter in 
the range < N,^ < 10 we obtain the constraints de- 
picted in Fig. ^ Joint contours in the c^jg — N,^ plane 
enclose the predicted value (1/3,3), and we do not find 
a significant degeneracy between the two parameters, 
since their impact on the CMB is almost orthogonal. 
As a consequence, the marginalized bound on c^jg is 
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FIG. 4: Relaxing the assumption on the background neu- 
trino density by including the number of neutrino fami- 
lies as a free parameter: joint 2-dimensional posterior 
probability contour plots in the cJig — Ni, plane, enclosing 
68% and 95% probability from the WMAP-hother CMB 
data alone (cyan/light gray), adding SLOAN matter power 
spectrum information (black, solid) and further imposing a 
conservative BBN constraint on N,^ (black, dashed). The 
indication that c^jg > is largely independent on prior 
knowledge of A^i^, since the joint WMAP-I-SDSS data can 
constrain c^jj, and N,y simultaneously. The prediction of 
the Standard Model is shown by the cross. 



only slightly weakened, and we obtain that c^jg = 
is disfavored at just above the 2a level (CMB and 
SLOAN data). Further imposing a conservative Big 
Bang Nucleosynthesis (BBN) constraint on N^, in the 
form of a Gaussian likelihood centered around = 3 
and with spread AN^, = 1 does not change the re- 
sult for c^jg in a significative way (see Fig. 

From a Bayesian model selection point of view, we 
can compare the evidence in favor of c^jg = 1/3 as op- 
posed to c^jg < 1/3 using the Savage-Dickey density 
ratio and taking a flat prior in the range < c^jg < 1/3 
(see '2^ for an explanation of the method and precise 
definitions). We obtain that WMAP-fSDSS data fa- 
vor Cyjg — 1/3 with odds slightly larger than 2:1, 
irrespective of the assumptions on the neutrino back- 
ground density, constituting positive (if only weak) 
evidence in favor of the Standard Model value. 

We also considered models with c^jg > 1/3, but in 
this case we found that the current data does not pro- 
vide any relevant constraint up to c^jg = 1. 

Conclusions In this letter we show that current 
cosmological data are providing for the first time an 
interesting indication for primordial anisotropics in 
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the cosmological neutrino background. Our result 
shows indication for the existence of a NB presenting 
anisotropics as predicted by the standard scenario, in 
which higher order multipoles in the neutrino distri- 
bution function are generated by free streaming. The 
significance is stih smah (shghtly larger than 2cr, quite 
independently on the assumptions on the radiation 
content) and possibly plagued by systematics. 

We analyzed a limited set of models and our con- 
clusions are valid only in the theoretical framework 
we consider. Enlarging the cosmological model by in- 
creasing the number of parameters could lower the 
significance of our result. For instance, a running of 
the spectral index would make the degeneracy be- 
tween cjjg and Us worse, thus weakening our con- 
straints. However, a measurable running is not ex- 
pected in the most common inflationary models and 
there is no strong indication for it from the data we 
considered, so it conservative to exclude this possibil- 
ity from our analysis. Another possibility would be to 
consider a dark energy model different from a cosmo- 
logical constant. The impact of dark energy is mainly 
in the position of the acoustic peaks and on the large 
scale CMB spectrum, which is rather orthogonal to 
the power suppression signature due to the viscosity 
parameter. Hence we do not expect a significant de- 
generacy between such models and c^jg. Finally, we 
made use of a very limited number of data-sets. In- 
clusion of Lyman-alpha data from SDSS, for example, 
would probably improve our constraints on the spec- 
tral index and therefore provide a tighter constrain 
on c^ig. The latest Ly-a analysis from SDSS 3] finds 



TT-s — 0.98 ± 0.02 which would further close our likeli- 
hood contours around the = 1/3 region. However, 
the precise statistical errors to be attached to such 
measurements are still under debate among special- 
ists, so we preferred not to include this data at the 
moment. 

One should note that the method we adopted pro- 
vides the only way to detect primordial anisotropics 
in the neutrino background in the foreseeable future. 
Furthermore, this work illustrates the power of current 
cosmological data-sets to test the hot Big Bang model 
and constrain subtle details of the energy density com- 
ponents of our universe. For instance, one might spec- 
ulate that detecting c^jg ^ 1/3 for the NB would hint 
to interaction of neutrinos with other species (one such 
example has been recently introduced by j2J| and in- 
vestigated in 113), including exotic components hkc 
dark energy. Future measurements and data releases 
should be able to strengthen our result, possibly shed- 
ding new light on the physics of the early universe. 
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